Tunneling between 2D electron layers with correlated disorder: anomalous sensitivity 

to spin-orbit coupling 
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Tunneling between two-dimensional electron layers with mutually correlated disorder potentials 
is studied theoretically. Due to this correlation, the diffusive eigenstates in different layers are 
almost orthogonal to each other. As a result, a peak in the tunnel I-V characteristics shifts towards 
small bias, V. If the correlation in disorder potentials is complete, the peak position and width 
are governed by the spin-orbit coupling in the layers; this coupling lifts the orthogonality of the 
eigenstates. Possibility to use inter-layer tunneling for experimental determination of weak intrinsic 
spin-orbit splitting of the Fermi surface is discussed. 

PACS numbers: 73.40Gk, 71.70.Ej, 72.25.Rb 
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I. INTRODUCTION 

Knowledge of spin-orbit (SO) splitting, A, of energy 
spectrum in 2D electronic systems is important for de- 
sign of spintronic devices in two respects. First, a num- 
ber of proposed schemes directly utilize the SO coupling 
for manipulating electron spin polarization by means of 
creating spatially inhomogeneous structures 1 . Second, in 
proposed schemes that are not based on SO splitting, the 
latter limits the device performance via a SO-induccd 
decoherence time 2 -^. Experimentally, large values of 
SO splitting can be extracted from conventional mea- 
surements, such as the beats of the Shubnikov-de Haas 
oscillations 4,5 . This, however, requires that At > 1, 
where r is the electron scattering time. Experimental 
determination of A in the opposite limit, At < 1, poses 
a considerable challenge. One has to look for physical ef- 
fects which are anomalously sensitive to the SO-coupling. 
An example of such effect is the weak localization/anti- 
localization crossover in magnetoresistance 6,7 . Tunnel- 
ing measurements offer another possibility. Even when 
At < 1, a structure related to A manifests itself in the 
I-V characteristics, provided that the disorder is long- 
range, so that AT tr > 1, where Tt r is the transport scat- 
tering time^. 

In 1993, Zheng and MacDonald 9 made an observation 
that, in the absence of the SO coupling, calculations of 
tunneling conductance between two parallel electron lay- 
ers with short-range but correlated disorder potentials is 
analogous to the calculation of conductance of a single 
layer with long-range disorder. Formally, both calcula- 
tions require solution of the equation for the vertex func- 
tions, obtained by a summation of ladder diagrams. For 
a single layer, the vertex function has a pole at frequency 
to = i/T tr , where 



of the vertex function is at u> = z/tq, where tq is defined 
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S(p 2 -\p+q\ 2 ), (1) 



where v = m/2n is the 2D density of states (per 
spin) and S(q) is the Fourier component of the corre- 
lator of the intralayer disorder potential, V(r); S(q) = 
J dr e _lqr (V (r)V (0)) . For interlayer tunneling, the pole 



— = 8ir 2 v 

T 



dq S(q)-S LR (q) 6(p 2 - |p + q| 2 ), (2) 



where similar to the above, Slr{q) is the Fourier compo- 
nent of the cross-correlator (Vl(t)Vr(0)) of the disorder 
potentials in the two layers. The physics captured by 
Eq. © is that despite strong scattering in each layer, 
the true eigenstates in both layers are almost identical 
when Vi, and Vr are strongly correlated. Then the pole 
at — iui — 1/tq <C 1/t reflects the fact that eigenstates 
in two layers with energy difference > 1/tq are almost 
orthogonal. 

Basing on the above analogy, pointed out by Zheng and 
MacDonald, one would anticipate anomalous sensitivity 
of the tunneling current between two layers with short- 
range correlated disorder to the SO splitting in the layers. 
In the present paper, we will illustrate this anomalous 
sensitivity for a particular example of tunneling between 
two identical quantum wells, to which electrons are sup- 
plied by a 5-layer of donors, located in the middle plane. 
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FIG. I: (A) Schematic illustration of a <5-layer of donors lo- 
cated symmetrically between two identical quantum wells. SO 
coupling constants in the wells have equal magnitude and op- 
posite signs; (B) Diagram describing tunneling current be- 
tween two wells with correlated disorder; (C) Diagrammatic 
equation for the vertex function. 
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II. TUNNELING CURRENT BETWEEN TWO 
2DEG WITH SPIN-ORBIT INTERACTION 

The system under study is shown in Fig. 1. Once 
the donors get ionized by yielding their electrons to the 
left and right electron gases, electric fields which they 
create in both layers are equal in magnitude and oppo- 
site in directions. As a result, the coupling constants 
in the SO Hamiltonians 10 of the two layers are opposite: 
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,(p x <x) z . The im- 



portant consequence of the geometry depicted in Fig. 1 
is that it allows to arrange correlation between spatial 
wave functions in different layers^ corresponding to dif- 
ferent energies separated by 2A. As a result, A manifests 
itself in the tunneling I-V characteristics. 
The tunneling Hamiltonian has the form, 



H 



'£ 



d 2 r ( v3i L)t (r)Vi K, (r) + ^(r)^>(r) 
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(3) 

where (r) and ^£ (r) are the electron operators in 
the two layers, and a is the spin index. The overlap inte- 
gral t for the size-quantization wavefunctions in the two 
layers is assumed to be real, for simplicity. The tunnel- 
ing described by Eq. ([3]) preserves both electron spin and 
momentum. 

Calculation of the interlayer tunneling current (see 
Fig. 1) reduces to finding the vertex function for the case 
when the electron Green's functions in the layers are ma- 
trices. Namely, the retarded Green's functions are 



6%\e, p) = [e - £ - a so (p x <r) 2 + i/2r}-\ 
G ( n\e, p) = [e-£ + a so (p x a) z + i/2r}- 1 



(4) 



where £ is the electron energy, measured from the Fermi 
level. Advanced Green's functions are obtained from 
Eqs. by reversing the sign of z/2r-terms. Solving the 
matrix equation illustrated in Fig. 1 yields the following 
generalized expression for the vertex 



T(w) = t 



{oj + i/rf 



4A 2 



(u + i/r)(u +i/r ) -4 A 2 



(5) 



where A = apF , and pf is the Fermi momentum. In the 
absence of the SO coupling Eq. © reduces to the result 
T(u)) = t(l — iojr) / (1 — ibJTo) of Ref. [9|- Incorporating the 
vertex function Eq. ((SJ) into the standard expression^ for 
the tunneling current 



I(V) 



{T(eV) 



e 2 AtV 3 



Tr / dp 6%\0,p)&f>{-eV,p)}, (6) 



we arrive to the final result, 

2e 2 t 2 AvV[4A 2 r- 1 + (e 2 V 2 + r" 2 )^ 1 ] 



I(V) = 



[e 2 V 2 - 4A 2 
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+ e 2 V 2 (r- 1 +r a 1 ) 2 



Here A is the lateral area. 

Anomalous sensitivity of the I-V characteristics ([7]) 
to the SO splitting is illustrated in Figs. 2-4. For large 
splitting, At > 1 (see Fig. 2), correlation between the 
disorder potentials is not important. The peaks of the I- 
V curves are located at eV = ±2A, while the peak widths 
are « 1/r. Such a form of the I-V curves reflects the fact 
that, with opposite signs the SO coupling constants in the 
layers, the intralayer spinor eigenfunctions are maximally 
correlated, when their energies differ by 2A. 




FIG. 2: The tunnel I-V characteristics is plotted from Eq. (5) 
for different values of dimensionless scattering rate At in the 
layers. For At > 0.5 the current is maximal at bias eV = 2A. 

As it is seen from Fig. 2, the position of the current 
maximum rapidly shifts from eV = 2A towards smaller 
biases for At < 0.5. The I-V characteristics for this case 
are shown in Fig. 3. A remarkable feature of the curves 
in Fig. 3 is their strong sensitivity to A when the disorder 
is strong, At < 1, so that the characteristics of individ- 
ual layers are insensitive to the SO coupling. For fully 
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At=0.15 
At=0.3 
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FIG. 3: The tunnel I-V characteristics is plotted from Eq. (5) 
for fully correlated disorder potentials in the layers at different 
values of At 1. The current is maximal at eV = 4A 2 t. 



(7) correlated disorder potentials in the layers, tq — > oo, the 
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I-V characteristics for different values of At < 0.2 fall on 
top of each other when plotted as a function of the ratio 
eV/AA 2 T. Thus, the position of maximum of the tunnel- 
ing current allows to extract the Dyakonov-Pereli^ spin 
decoherence time t s — (2A 2 t) _1 . The underlying reason 
is that, due to opposite signs of the intralayer SO cou- 
pling constants, the eigenfunctions in the layers are not 
orthogonal even if disorders are fully correlated. Then 
T.T 1 is a quantitative measure of the energy interval in 
which the orthogonality is lifted. The fact that position 
of the maximum in Fig. (3) is at eV = 2/r s reflects that 
electrons in both layers undergo spin relaxation. 

Incomplete correlation of disorder potentials, Vt,(r) 
and Vu(r), in the layers is another source of lifting of 
orthogonality of eigenstates. This mechanism is quanti- 
fied by the energy scale 1/to, defined by Eq. ^j. It might 
be expected that in the presence of both mechanisms the 
maximum is located at eV = 2/r s + 1/to, which has a 
meaning of a combined dephasing time. This is indeed 
the case, as illustrated in Fig. 4. 



III. EFFECT OF ELECTRON-ELECTRON 
INTERACTIONS 

Let us address the question whether the above SO- 
induced peaks survive the presence of electron-electron 
interactions. Interactions cause a dynamic lifting of or- 
thogonality of the eigenstates, and might result in the 
broadening of the peaks. We now demonstrate that at 
zero temperature the peaks are robust, but eventually 
get smeared away as the temperature increases. 
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FIG. 4: I-V curves for different values of spin-orbit coupling 
and different degrees of correlation of disorder in the layers. 
The maxima are at biases determined by the combined de- 
coherence rate. The latter is dominated by 1/ro term for 
At = 0.05 and by spin relaxation term, 4A 2 r, for At — 0.4. 

On the quantitative level, in order to incorporate both 
the interactions and the correlated disorder into the the- 
ory, it is convenient to express the tunneling current in 
terms of the exact eigenfunctions, which are the same in 



the two layers. Let us denote with ip m ( r ) the Tn-th eigen- 
state for a given realization of disorder potential. The en- 
ergy of this state is equal to e m as electron-electron inter- 
actions are neglected. In the presence of electron-electron 
interactions the retarded electron Green function can be 
written as, 



G R (e, ri,r 2 ) = 



jMri)V4(r2) 
e — £m — S m (e) 



(8) 



where S m (e) denotes the electron self-energy of the m-th 
eigenstate. 

The knowledge of the eigenfunctions suffices to evalu- 
ate the tunneling current (in the lowest order in t) in a 
general form, 



I{V) = lei 2 



2vr 



n(e) — n(e + eV) 



J dx t dx 2 A^ (e, r ls r 2 )AW ( e + eV, r 2 , n), 



(9) 



and express it via the Fermi-Dirac distribution n(e) 
and the non- averaged spectral functions in the left and 
right layers, A- L \e, T\, r 2 ) and A^ R > (e, ri, r 2 ), respec- 
tively. The spectral function is determined by the dif- 
ference of the retarded and advanced functions in each 
layer, 



A(e,r 1 ,T 2 ) = - \G R (e,ri,r 2 ) 



G A {e,ri,r 2 ) 



(10) 



At this point, we emphasize that, for fully correlated dis- 
order (neglecting spin-orbit interaction), we can perform 
coordinate integration in Eq. ^ prior to performing the 
configuration averaging, and cast it in the form 



u 

I(V)=2et^ J 



dc 
2^ 



1 



-eV 



6 e m S m (e) 
1 



e + eV - e m - £ m (e + eV) 



(11) 



The reason why explicit integrations over ri and r 2 can 
be performed in Eq. ([9]), leading to Eq. (fTTj) . is the 
mutual orthogonality of the eigenstates in the two layers 
resulting from the fully correlated disorder. In Eq. (fTT|) 
we have also set T = in the difference of the Fermi 
functions. This is justified as long as T is much smaller 
than Ep. The reason is that, similarly to the in-plane 
conduction, the temperature dependence of the tunnel 
current comes exclusively from the T-dependence of the 
self-energy, i.e. from inelastic processes. 

The question whether £ m (e) and S m (e + eV) in 
Eq. (fTTj) can be replaced by the disorder averaged val- 
ues is highly non-trivial in the limit T — > 0. However, 
we can rigorously address the issue of smearing of the 
SO-related peak in the I-V curve for disordered layer by 
treating interactions at the perturbative level, which cor- 
responds to the expansion of Eq. (fTT|) to the first order 
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in £ m . This expansion yields, 



6I(V) 



2t 2 vA 
eV 2 



de 



lc {e + eV)+ lt+eV (e) , (12) 



where 7 e (w) = v 1 X) m (^( e ~ e m) SS m (cj)) is now the 
disorder-averaged inverse inelastic lifetime. 

(i) For At <C 1 the peak position, eV — 4A 2 r, is below 
the elastic scattering rate, i.e. at the energy correspond- 
ing to the peak position the motion of electrons is diffu- 
sive. The corresponding lifetime was studied in the sem- 
inal paper Ref.[l4|, and was shown to be, 7(e) ~ \e\/E F T. 
We can now compare Eq. (TT2"]) with the "non-interacting" 
value of current, given by Eq. ([7]), at the bias correspond- 
ing to the peak position, eV — 4A 2 r. We find that the 
ratio 81 /I ~ 1/ E f t is small regardless of the actual value 
of the decoherence rate. 

(ii) Similarly, for large values of spin-orbit splitting, 
At 3> 1, we should utilize ballistic inverse lifetime, 
7(e) = (e 2 /4wE F )ln[E F /e\, established in Refe. USES 
Comparison of Eq. (TT2"|) with the value given by Eq. (J7J at 
the peak position, eV — 2A, we conclude that the corre- 
sponding ratio is again small, SI /I ~ (1/E f t) In [E F /A]. 

This suggests that, at zero temperature, interactions 
do not destroy the SO-induced peak in the I-V curve. 
However, this destruction eventually happens upon in- 
creasing T. A crude estimate for the temperature at 
which the peak is washed out by interactions can be 
obtained by equating the peak position eV = 4A 2 t to 
7(e = T). With logarithmic accuracy, this yields the re- 
striction T < (At) 2 E F: so that even with At < 1 the 
peak survives at reasonably high temperatures. 

To trace quantitatively the smearing of the peak with 
T, we first note that "single-electron" I-V characteris- 
tics (J7J can be obtained from Eq. (fTTj) upon inserting 
spin decoherence rate into the self-energy, E TO — * 2iA 2 r, 
and replacing the sum over e m — > £ by the integral, 
Em — * v ^ I As the next step, we take into account 
the fmite-T decoherence by writing £ = 2iA 2 r + i-fr, 
where 7T = (T/2E f t) In (Ti/T) 14 - 17 and T x = r 2 E F T 3 ; 
here r s is the interaction parameter of 2D electron gas. 
To utilize the energy-independent j(T) in the self-energy, 
the temperature must be large compared to the peak po- 
sition, eV — 2A 2 t. This requirement does not contradict 
the restriction on the smearing obtained from the above 
crude estimate. Indeed, both conditions can be conve- 
niently rewritten as, T/E f t <C A 2 r <C T, so that it is the 
large value of E f t which makes them consistent. Upon 
the suggested replacements, the temperature-dependent 
I-V characteristics follows from Eq. (J7J with the spin 
relaxation rate modifies as 2A 2 r — > 2A 2 r + 77-, with 

T = OO 



I(V) = Ae 2 t 2 Av 



V(2A 2 t + 1t ) 
e 2 U 2 + 4(2A 2 T + 7T ) 



(13) 



Eq. (TT3"|) indicates that the position of maximum of the 
I-V curve shifts almost linearly with temperature, see 



Fig. 5. This suggests that the SO relaxation rate can 
be inferred from experiment even when measurements 
are performed at T > (At) 2 E F - One has to plot the 
peak position as a function of T and extrapolate the 
data to T — > 0. Besides, the actual restriction on T is 
"softer" than the one obtained from the crude estimate, 
by virtue of numerical coefficients in 77^ and spin relax- 
ation time 1/t s . Indeed, the requirement 7tt s /2 < 1 
imposes (neglecting logarithmic factor) the following re- 
striction, T < A(At) 2 E f . 



I(V) / (4et z Av) 



At=0.05 
At=0.2 




FIG. 5: (Color online) Temperature dependence of the tun- 
neling characteristics for different values of spin-orbit coupling 
constant and fully correlated disorder, from Eq. (|13p . Increas- 
ing T results in the broadening of the peak and its net shift 
towards higher biases. 



IV. SUMMARY AND CONCLUSIONS 

Our main finding is that, with correlated disorder in 
the layers, the SO coupling causes a zero in dl/dV even 
for At -c 1, when the spin subbands in the layers are 
not resolved. For clean layers with At > 1, sensitivity of 
tunneling current to the SO coupling was pointed out in 
Ref. [H 

The condition that position of zero in dl/dV is due 
to the SO coupling is that the contribution, 4A 2 t, to 
the combined decoherence rate exceeds 1/to, caused by 
incomplete correlation of disorders in the layers. To 
estimate the feasibility to meet this condition, we as- 
sume that the origin of incomplete correlation is a fi- 
nite width, a, of the <5-layer, see Fig. 1. Assuming 
that the in-plane positions of donors with concentra- 
tion, Nd, are completely random, the Fourier compo- 
nents of the correlators, entering the expression Eq. @ 
for Tq 1 , can be presented as Slr.(h) = A r d|£/(q)| 2 e -9 ' i ; 
5(q) = A^|£/(q)| 2 sinh (qa) e~ qd /qa, where d is the bar- 
rier thickness, and U(q) is the Fourier transform of the 
potential created by a donor in the layer. Then Eq. @ 
takes the form 



— = 2ttza\, 
to 



d J dq \U(q)\ : 



( sinhqa 

V qa 



- 1 



(14) 
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Assuming that the screening radius in the layers is 
smaller than d, we can set U(l/d) « U(0). Then Eq. (TT4"|) 
yields ro/r = (a/d) 2 , so that the condition 4A 2 t > Tq -1 
reduces to At > a/2d, i.e. the value At = 0.05, used in 
the numerics above, is quite feasible for the atomically 
sharp 6- layer. 

As a final remark, the assumption, crucial for our cal- 
culations, was that the barrier is spatially homogeneous, 
so that the tunneling occurs with the conservation of the 
in-plane momentum. We had also assumed that the posi- 
tions of donors in the (5-layer are random. This random- 
ness might, in principle, lift the momentum conservation. 
The condition that the effect of randomness is negligible 
is that the tunneling-induced splitting, t, of the states 
in the layers is the smallest scale in the problem. In 
fact, we had used this condition by restricting the cal- 



culations to the lowest order in t. Under this condition, 
the under-barrier scattering of the tunneling electron by 
donors in the classically forbidden region is exponentially 
suppressed as compared to the situation when donors are 
located in the vicinity of the layers. 
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